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Abstract. A detailed spectroscopic study of the ISO- 
SWS data of the red giant a Tau is presented, which en- 
ables not only the accurate determination of the stellar 
parameters of a Tau, but also serves as a critical review 
of the ISO-SWS cahbration. 

This study is situated in a broader context of an itera- 
tive process in which both accurate observations of stellar 
templates and cool star atmosphere models are involved 
to improve the ISO-SWS calibration process as well as the 
theoretical modelling of stellar atmospheres. Therefore a 
sample of cool stars, covering the whole AO - M8 spectral 
classification, has been observed in order to disentangle 
calibration problems and problems in generating the the- 
oretical models and corresponding synthetic spectrum. 

By using stellar parameters found in the literature 
large discrepancies were seen between the ISO-SWS data 
and the generated synthetic spectrum of a Tau. A study 
of the influence of various stellar parameters on the theo- 
retical models and synthetic spectra, in conjunction with 
the Kolmogorov-Smirnov test to evaluate objectively the 
goodness-of-fit, enables ns to pin down the stellar pa- 
rameters with a high accuracy: Tcfr = 3850 ± 70 K, log 
g= 1.50±0.15, M 2.3±0.8 M©, z = -0.15±0.20 dcx, 
= l.TibO.a kms-i, = 10±1, e(C) = 8.35±0.20 

dex, e(N) = 8.35±0.25 dex, e{0) = 8.83±0.15 dex and 0d 
= 20.77 ± 0.83 mas. These atmospheric parameters were 
then compared with the results provided by other authors 
using other methods and/or spectra. 
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1. Introduction 

The modelling and interpretation of the ISO-SWS (In- 
frared Space Observatory - Short Wavelength Spectrome- 
ter) data require an accurate calibration of the spectrom- 
eters (Schaeidt et al. 1996). In the SWS spectral region 
(2.38 - 45.2 fim) the primary standard calibration candles 
are bright, mostly cool, stars. The better the behaviour 
of these calibration sources in the infrared is known, the 
more accurate the spectrometers can be calibrated. ISO 
offered the first opportunity to obtain continuous mid- 
infrared spectra between 2.38 and 45.2 fim at a spectral 
resolving power of ~ 1500, not polluted by any molecu- 
lar absorptions of the earth's atmosphere. Therefore our 
knowledge on the mid-infrared behaviour of the stellar 
calibration sources is limited. Refining the synthetic refer- 
ence spectra used to calibrate the SWS can only be done 
by refining the model atmospheres of the stars. A full ex- 
ploitation of the ISO data can therefore only result from 
an iterative process in which both accurate observations 
and new modelling are involved. In order to obtain a re- 
liable convergence in such an iterative process, one needs: 
(1) a sample of bright stellar objects spread over a large 
range in spectral type; (2) a thorough understanding of 
cool star atmospheres and of the influence of the vari- 
ous parameters — e.g. Toft, log g, and chemical composi- 
tion — on the emergent spectra. With this in mind, stars 
with spectral types ranging from AO - M8 were observed 
with ISO-SWS. It is important to cover a broad parame- 
ter space in order to distinguish calibration problems from 
problems related to the model and/or in the generation of 
the synthetic spectrum. Concerning point 2 above, it is 
necessary to perform an intensive study of the influence 
of different stellar parameters on synthetic spectra. The 
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results of such a study will be presented with the red gi- 
ant a Tau as a test case. The results for 15 other stars in 
our sample will be presented in forthcoming papers. 

This research has been done within the framework of 
the proposal STARMOD (Title: Accurate modeUing of 
cool-star atmospheres; P.I.: C. Waelkens; C.I.: M. Cohen, 
L. Decin, Th. de Graauw, L.B.F.M. Waters) and the ZZ- 
STARM proposal (Title: Accurate modelling of cool-star 
atmospheres; P.I.: L. Decin; C.I.: M. Cohen, C. Waelkens, 
Th. de Graauw, L.B.F.M. Waters). Some calibration data 
have been provided by the SWS Instrument Dedicated 
Team (SIDT) in the framework of a quick-feedback refin- 
ing of the model Spectral Energy Distribution (SED) of 
the calibration sources used for the SWS calibration. 

So far, the analysis of the discrepancies between the 
ISO-SWS data and the corresponding synthetic spectra 
has been restricted to the wavelength region from 2.38 
to 12 ^m, since the lack of comprehensive molecular and 
atomic line lists hamper fast progress at longer wave- 
lengths (12 - 45 [iva). Furthermore the brightness of the 
stars drops quickly in this wavelength region so that the 
same signal to noise ratio will not be achieved. A third 
point is that the SEDs may also be affected by unknown 
circumstellar contribution. 

This paper is organized as follows: in Sect. 2 the obser- 
vations are described and the data reduction procedures 
are discussed. In Sect. 3 a summary of the literature con- 
cerning a Tau is presented, on the basis of which the 
starting values for our modelling stellar parameters are 
selected. The effect of changing stellar parameters on the 
synthetic spectra of K and M giants is analyzed in Sect. 
4. In Sect. 5 the method of analysis is outlined, while in 
Sect. 6 the results are discussed. In the last section, Sect. 
7, the main conclusions are summarized. 

2. Observations and data reduction 

A full scan observation of a Tau has been obtained in 
the 2.38 - 45.2 ^m wavelength range with the Short- 
Wavelength Spectrometer (SWS) (de Graauw et al. 1996) 
on board ISO. This has been done by use of the SWS ob- 
serving mode AOTOl (= a single up-down scan for each 
aperture/order combination) with scanner speed 4, result- 
ing in a resolving power of ~ 1500. The observation lasted 
for 6650 s and was performed during revolution 636. 

The data were processed to a calibrated spectrum sam- 
pled in all the pixels, the so-called Auto- Analysis product, 
using the procedures and calibration files of the ISO off- 
line pipeline version 7.0. For a description of the flux and 
wavelength calibration, we refer to Schaeidt et al. (1996) 
and Valentijn et al. (1996). This reduction resulted in an 
oversampled spectrum for each of the four AOT-bands, 
which each consist of 12 detectors. 

Remaining instrumental effects, such as fringes in the 
12.0 - 29.5 /im part of the spectrum, were removed using 
the SWS Interactive Analysis package (lA) provided by 



the SIDT. The band-2 data (4.08 - 12 ^m) are severely 
affected by detector memory effects. Time and flux depen- 
dent memory effects occur for all four instruments (SWS, 
LWS, PHOT, CAM) on board ISO. So far, several at- 
tempts have been made to model the corrections of the 
transient effects (see e.g. Van Malderen et al. 1999), but 
currently no foolproof method exists to correct the SWS 
data for these transient effects. Memory effects appear to 
be less severe in the down-scan measurements, which are 
obtained after the up-scan data, suggesting a more stabi- 
lized response to the flux level for down-scan data. The 
SWS flux calibration of band 2 was corrected for a regime 
where the flux history of the detectors is comparable to the 
second scan in the observation, the so-called down scan. 
We therefore used the down-scan data of our observation 
as a reference to do a correction of the flux level of the 
first scan. 

Several procedures were performed on each sub-band 
separately to combine the twelve detectors. A sigma-clip- 
ping procedure (with tr = 2.0) was used to reduce the noise 
by discarding the datapoints in every resolution element 
for which the difference to the mean was > 2(t. After align- 
ing the twelve different detector signals to their average 
level and rebinning to the expected resolution (see Table 
1), the final result was obtained. Although the resolution 
changes in within one sub-band, the resolution was kept 
constant (see Table 1) and the value was taken to be the 
lowest value in Fig. 6 by Lorente (1998). The individual 
sub-band spectra, when combined into a single spectrum, 
can show jumps in flux levels at the band edges. This is 
due to imperfect flux calibration or wrong dark-current 
subtraction for low-flux observations. Using the overlap 
regions of the different sub-bands and looking at other 
SWS observations, several sub-bands were multiplied by a 
small factor (see Table 1) to construct a smooth spectrum. 
Due to the problems with memory effects in band 2, the 
factor of band 2b was determined by use of the template 
of Cohen et al. (1992) and the obtained synthetic spectra. 
Note that all shifts are well within the photometric ab- 
solute calibration uncertainties claimed by Schaeidt et al. 
(1996). 

3. Literature study: a. Tau 

a Tau (= HD 29139 = HR 1457 = HIC 21421) is classiflcd 
as a K5III giant. A detailed literature study was necessary, 
on the one hand to extract the best possible set of starting 
fundamental parameters in order to reduce the number of 
calculated spectra and on the other hand to check the con- 
sistency between the stellar parameters deduced from the 
ISO-SWS data and other spectra/methods. An exhaustive 
discussion of published results is presented in Table 2 and 
in Appendix A. 

After evaluation of the various sources, the stellar pa- 
rameters of McWflham (1990), with Tes = 3900 K, log g 
= 1.60, it = 2.0 kms~\ [Fc/H] = -0.30, together with 



Table 2. Literature study of q Tau, with the effective temperature Tofi given in K, the mass M in Mq, the microturbulent velocity in kms^ , the angular diameter 
Od in mas, the luminosity L in Lq and the radius R in Rq. Angular diameters deduced from direct measurements (e.g from interferometry) are written in italic, while 
others (e.g. from spectrophotometric comparisons) are written upright. 
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Table 1. Resolution and factors used to shift the sub-bands. 
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e(C) = 8.27, e(N) = 7.86, e(0) = 8.60 (Lambert & Rics 
(1981), but with [N/Fe] taken to be -f0.20 according to 
Luck & Challener (1995), see Appendix A), a mass of 1.5 
M0 (van Paradijs & Meurs 1974), the limb-darkened an- 
gular diameter 9d = 20.88 mas from Ridgway et al. (1982) 
and = 10 (Smith & Lambert 1985) were taken 

as starting values. In Fig. 1 both band 1 of the ISO-SWS 
data of a Tau and the synthetic spectrum generated with 
these parameters are shown. From Fig. 1, it is already 
obvious that not only the angular diameter 6d should be 
lowered and/or Toft should be increased, but also other 
parameters may need improvement. 



Alpha Tau 




2,5 3.0 3 5 .1 

Wovelength [n-icron] 



Fig. 1. Comparison between band 1 of the ISO-SWS data of 
a Tau (black) and the synthetic spectrum (grey) with stellar 
parameters T^s = 3900 K, log g = 1.60, M = 1.5 Mq, z = 
-0.30, = 2.0 kms~\ ^^C/^C = 10, e(C) = 8.27, e(N) = 
7.86, £(0) = 8.60 and Oa = 20.88 mas. 



4. Influence of stellar parameters on synthetic 
spectra 

In this section the effect of stellar parameters on the ab- 
sorption by CO, SiO, OH and H2O is discussed. We mainly 
focus on these molecules since they are the most promi- 
nent absorbers for oxygen-rich giants in the wavelength 
range considered here (2.38 - 12 fj.m) (see, e.g., Fig. 2 and 
Fig. 3 in Decin et al. (1997) and Fig. 4 in this paper). 
The goal of this part of the study is to learn how the syn- 
thetic spectrum will change when one of the parameters 
is changed within its uncertainties, the only exception be- 
ing the stellar mass and the ^^C/^^C ratio due to their 
small influence. A simple equation for any chemical com- 
pound is not easily obtained, nor is a unique scenario able 
to explain all different behaviours, but some patterns do 
emerge. These results are then useful for the determination 
of the origin of the discrepancies seen between ISO-SWS 
data and synthetic spectra. 

Since the strength of molecular (and atomic) lines is 
proportional to the ratio of line to continuous absorption 
coefficient, approximations for this ratio are sought 

for. The approximate ratio is referred to as TZ and equa- 
tions for some molecules have been deduced according to 
the approach of Kjaegaard et al. (1982) (see Appendix B). 

The ratio TZ has been studied for the surface flux 
using the following parameters: 
Tefi' = 3650 K - 3850 K - 4050 K 
log g = 1.00 - 1.50 - 2.00 
z = -0.15 - 0.00 - 0.15 
M = 1.5 Mq - 10 Mq - 15 Mq 

= 1.0 kms-i - 2.0 kms"^ 
12C/13C = 10 - 4.26 
£(C) = 8.24 - 8.54 
£(N) = 8.26 - 8.56 
£(0) = 8.83 - 9.13 

All computed spectra are compared with the synthetic 
spectrum with parameters: Teff = 3650 K, log g = 1.00, 
z = 0.00, M = 1.5 Mq, = 2.0 kms-\ ^^C/^^C = 10, 
£(C) = 8.24, £(N) = 8.26, £(0) = 8.83. If not specified 
these stellar parameters are used. In Fig. 2 the overall 
results of a change ATeff = 200 K (a), Alog g = 0.50 (b), 
AM = 13.5 Mq (c), Az = 0.15 (d), A^t = 1.0 kms"! (e), 
A12C/13C 5.84 (f), A£(C) = 0.30 (g), A£(N) = 0.30 (h), 
A£(0) — 0.30 (i) are shown. The results mentioned and 
discussed in this study apply only to models with similar 
stellar parameters. 

4.1. Models 

The models and corresponding synthetic spectra have 
been computed by using the MARCS-code (Gustafsson et 
al. 1975). Since 1975, this code has undergone some mod- 
ifications, the most important ones being the replacement 
of the Opacity Distribution Function (ODF) technique by 



L. Decin et al.: ISO-SWS and modelling of cool stars 



2,4C 2,45 2,50 2,55 2-60 2,7C 2,80 2,90 5,00 5,'.C 5,20 5,50 5,40 5,50 5,60 3,70 5,80 5,90 4,00 



4,5 4,8 5,0 5,2 5,5 6,0 



7,0 3,0 9,0 10,0 11,0 12,0 



1,00 

0,95 



1,00 
0,95 



1,05 
1,00 
0,95 

?:?8 

1,05 



1,00 

0,95 



1,00 
0,95 









SIO 
kv=2 

\ 

T^,=3650 K / T,„=3850 K ' 


CO 
/ 




(«) 

Sic 








CO 
Ai'=2 


HjC + CH 




SiO 
Ai'=2 ■ 


CO 

\ 


H;0 + OH 


(l>) 

SiC 

ai,=i 




log g-1.0 / log g-1.5 




















■ — — 


— 






U=1.5 Me / M=15 














2=0.00 / z=0.15 














CO 
4v=2 






SiO 
Av=2 

'f' 

J,-2.0km/s/ f,-l,Okm/s 


CO 

Av=l 




(=) 

SiO 


CO 
Av=2 












!") 






'*C/'^C-10.1 / 'V*C-4.26 


CO 
Sv=l 






CO ^ 
Av=2 












(g) 




£(C}=8.24 / e(C)=8,54 








CN 




















£(N)=S-26 / e(N)=8.55 










HjO + OH 














e(0}-8.83 / e(0)-9.13 









2.4C 2,45 2,50 2,55 2,50 2,70 2,80 2,90 3,00 3,10 3,20 3,30 3,40 3,50 3,60 3,70 3,S0 3,90 4,00 

wavelerglh [fim] 



4,6 4,B 5,0 5,2 5,5 



7,0 8,0 9,0 10,0 



Fig. 2. Division of tlie syntlietic spectrum witli parameters Tcff = 3650 K, log g = 1,00, M = 1,5 M0, z = 0,00, — 2.0 kms ^, 
12q^13q = 10, e(C) = 8.24, e(N) = 8,26, e(0) = 8,83 by the synthetic spectrum with the same parameters but with (a) T^ff 
= 3850 K, (b) log g = 0.50, (c) M = 15 M©, (d) z = 0,15, (e) = 1 kms-\ (f) i^c/^C = 4,26, (g) e(C) = 8,54, (h) 
e(N) = 8,56 and (i) e(0) = 9,13, Prominent molecular bands are specified in the plot. Plot (a) has been divided by the factor 
3650/3850 to compensate for the higher flux of the spectrum with an effective temperature of 3850 K, 



the Opacity Sampling (OS) technique, the possibihty to 
use a spherically symmetric geometry for extended objects 
and major improvements of the hne and continuous opac- 
ities (Plez et al. 1992). 

The common assumptions of spherical or plane-parallel 
stratification in homogeneous stationary layers, hydrostat- 
ic equilibrium and LTE were made. Energy conservation 
was required for radiative and convective flux, where the 
energy transport due to convection was treated through 
a local mixing-length theory (Henyey et al. 1965). The 
mixing-length I was chosen as 1.5 iJp, with Hp the pres- 
sure scale height, which is a reasonable quantity to sim- 
ulate the temperature structure beneath the photosphere 



(Nordlund & Dravins 1990). Turbulent pressure was ne- 
glected. The reliability of these assumptions is discussed 
by Plez et al. (1992), 

The synthetic spectra were generated using the Tur- 
boSpectrum program described by Plez et al, (1993), and 
further updated. The program treats the chemical equi- 
librium for hundreds of molecules with a consistent set of 
partition functions and dissociation energies. Solar abun- 
dances from Anders & Grevesse (1989) have been as- 
sumed, except for the iron abundance, e(Fe)= 7.51, which 
is in better agreement with the meteoritic value. 

The continuous opacity sources considered here are 
H", H, Fe, (H+H), H+, H^, He I, He Iff, He", C I, C H, 
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C Iff, C Ilff, C-, N I, N II, N-, O I, O II, 0", CO", 
H2O-, Mg I, Mg II, Al I, Al II, Si I, Si II, Ca I, Ca II, 
H2(pr), Ile(pr), e^^, H^^, Il2sc, where 'pr' stands for 'pres- 
sure induced' and 'sc' for 'scattering'. 

For the line opacity in the SWS range (2.38 - 12 /im) a 
database of infrared lines including atoms and molecules 
has been prepared. For the atomic lines the data listed 
by Hirata & Horaguchi (1995) have been used, for CO 
those by Goorvitch (1994), for SiO those by Langhoff & 
Bauschlicher (1993), for II2O those by J0rgcnscn (1994) 
and Ames (Partridge & Schwenkc 1997), OH lines by Sau- 
val (Melen et al. 1995), Schwenkc (1997) and Goldman ct 
al. (1998), NH lines by Sauval (Grevesse et al. 1990; Geller 
et al. 1991) and CH lines by Sauval (Melen et al. 1989; 
Grevesse et al. 1991) and CN lines by Plez (private com- 
munication). The dissociation energy for CN was taken 
to be 7.76 eV. An exhaustive discussion on the accuracy 
and completeness of infrared spectroscopic line lists can 
be found in Dccin (2000). From this study, a preference 
emerged for the II2O line list by Ames and the OH line 
list by Goldman. To remove the small difference between 
the originally 'vacuum' wavelengths for the ISO observa- 
tions and 'air' wavelengths in spectroscopic linclists above 
200 nm, Edlen's formula (Edlen 1966) was extended to the 
infrared. 



Therefore a given line (J) intensity could be increased or 
decreased according to its J-value only. 




-2 2 A -8 -b -4 -2 



Fig. 3. Comparison between the partial pressure of several 
molecules for the models with TcH — 3650 K (full line) and 
Toff = 3850 K (dashed line) and with log g = 1.00, M = 1.5 
Mq, z = 0.00, C* = 2.0 kms-\ ^^c/i^c = 10, e(C) = 8.24, 
e(N) = 8.26, e(0) = 8.83. In plot (a) the partial pressure of 
the molecules is given, while in plot (b) the partial pressure of 
the molecules is divided by p(H I)*Pc. 



4-2. Effect of changing the effective temperature 

When the temperature increases, fewer molecules are for- 
med resulting in a smaller line absorption coefficient. At 
2.3 iim, Kv is primarily due to H~ free-free absorption, 
thus its value depends directly on the electron pressure Pe- 
Increasing the temperature causes more ionization events 
and thus more free electrons, but at the same time the H~ 
ion itself is less easily formed. Fig. 3a. shows that, for the 
models under consideration, this latter effect is the most 
important one and so decreases when the effective tem- 
perature is increased from 3650 K to 3850 K. Moving in- 
wards from the outer photosphere, the line absorption co- 
efficient li/ of all molecules (not H~ ) reaches its maximum 
at the location where the effect of an increasing temper- 
ature overtakes the effect of a higher density. The partial 
pressure p of H~, being proportional to p(H I)*Pe, keeps 
on rising because of the increase of the number of free elec- 
trons at higher temperature. When the line-to-continuum 
contrast has to be known, one has to consider the partial 
pressure of the molecule divided by (p(H I)*Pe) (see Fig. 2 
and Fig. 3b). The simple statement that there are weaker 
absorption bands for higher temperatures does not always 
hold: the approximate ratio TZ of CO, SiO, OH and H2O 
decreases, but for CN the lower continuous opacity com- 
pensates for the lower line opacity (see Fig. 3b). Further- 
more, the relative population of rotational levels strongly 
depends on T (the maximum population occurs at a Jmax- 
value and a lower T corresponds to a lower Jmax- value). 



4.3. Effect of changing the gravity and the mass 

Although g and M are closely related, the effect of chang- 
ing the mass is smaller than that of changing the gravity, 
because the first one only affects the extension and the 
latter one also changes the pressure structure of the atmo- 
sphere. This is clearly visible in Fig. 2. The assumption of 
hydrostatic equilibrium links the gas pressure to the sur- 
face gravity. In a cool star, the gas and electron pressure 
(Pg and Pe) can be written as: 

Pg « constant ■ gP (0.59 <p< 1.05) (1) 
Pe sa constant ■ gP (0.20 <p< 1.03). (2) 

Due to the higher extension of the atmospheres considered 
here, this value of p has a larger range compared to the 
value discussed by Gray (1992). 

When the gravity increases, the electron pressure in- 
creases rapidly (Eq. 2), resulting in a higher H~ free- 
free absorption. Other effects of increasing the gravity are 
higher number densities supporting the molecular forma- 
tion and the fact that molecules can be formed till higher 
temperatures, both increasing l^. The final result is that 
the approximate absorption coefficient TZ can either in- 
crease or decrease depending on the relative change of l^, 
and K^. 

Using the approximate relations according to Kjasr- 
gaard et al. (1982) (see Appendix B), the conclusion is 
reached that an increase of gravity leads to a decrease 
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of the strength of the CO, NH and CN lines and to an 
increase of the strength of the H2O hnes, with no depen- 
dence for the OH and SiO lines. The increasing gravity 
does lead to a strengthening of the H2O lines and also of 
the OH lines while it diminishes the 1^1 Hi> ratio of CO, 
NH, SiO and CN (Fig. 2). From the approximate rela- 
tions for TZ no quantitative conclusion on the parameter 
dependence can thus be made. 

An increase in mass yields the inverse (but smaller) 
effect on the partial pressure (and so on the line absorption 
coefficient l^) as an increase in gravity. 

4.4. Effect of changing the metallicity 

For the computations with a higher metallicity, [atom/H] 
has been increased from 0.00 dcx to 0.15 dex for all atoms 
except for H, He, C, N and O. Because an increase in 
metal abundance leads to a decrease in gas pressure at 
each optical depth (as k^, increases), the partial pressures 
of the molecules CO, OH, N2, H2O and CN consequently 
decrease, while p(SiO) and p(TiO) are somewhat higher 
because of the higher abundance of Si and Ti respectively. 
An increase in the overall metal abundance increases the 
amount of electrons (Pe) and as a consequence also the 
value of Ki, formed by H~ . This increase of k^, dominates 
the changes in the spectrum: all molecular features be- 
come fainter at increased metallicities (with C, N and O 
kept constant) (Fig. 2) . This can also be seen in the equa- 
tions for TZ (see appendix B) for the extreme case in which 
all carbon is present in the form of CO with Toff < 4500 
K. The strong metallicity dependence of H2O in the ap- 
proximate relation for TZ may be well visible, though this 
increase (with decreasing z) is small for weak H2O bands. 

4-5. Effect of changing the microturhulent velocity 

The impact of changing the microturhulent velocity is 
most important for the CO lines with a large equivalent 
width. When a line is saturated, increasing widens the 
wavelength range covered by the absorption and reduces 
the saturation, thus increasing the total absorption. In the 
line center of an unsaturated line, a smaller microturbu- 
lence corresponds to a higher absorption coefficient at that 
frequency (see Fig. 2). 

Decreasing the microturhulent velocity from 2 kms^^ 
to 1 kms~^ diminishes the CO first overtone absorption by 
maximum 10% and the CO fundamental by up to 6 %, the 
SiO first overtone by up to 4%, and the SiO fundamental 
by up to 3%, the fundamental band of OH by up to 3% 
and the ro-vibrational stretching modes {vi and v^) and 
bending mode (1^2) of H2O by up to 0.5%. 

4-6. Effect of changing the isotopic ratio ^^C/^'^C 

Red giant atmospheres display altered C, N, and O abun- 
dances and isotopic ratios (in particular ^^C/^'^C), due 



to dredge-up episodes. Suppose that the ^^C/^'^C ratio is 
changed from 10 (case I) to 4.26 (case II). If a ^■^CO line is 
saturated in case I, its equivalent width will decrease very 
few from case I to case II. If we assume the ^'^CO line to be 
very weak in case I, its equivalent width will increase very 
much from case I to case II by a factor nearly equal to the 
new/old isotopic ratio. Therefore it is understandable that 
a decrease in ^^C/^^C will consequently cause 7^(C0) to 
increase, although this increase may be very small if the 
CO bands are weak (Fig. 2). 

4-7. Effect of changing e(C), e(N) or e(0) 

When the ratio C/0 increases (but remains < 1), the 
structure of the atmosphere is affected through changes 
in line blanketing. More CO will be formed and less oxy- 
gen will be left (after CO formation) to form oxygen-based 
molecules. Hence, as illustrated in Fig. 2, the total absorp- 
tion of CO will increase, while it will decrease for SiO, OH 
and H2O. Increasing the oxygen abundance gives roughly 
the reverse effect, but the total CO absorption hardly 
changes. The largest effect of an increasing nitrogen abun- 
dance is an increase of the NH and CN line strengths. 

5. Method of analysis 
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Fig. 4. Relative importance of atoms and molecules in a model 
with Teff = 3850 K, log g = 1.00, M = 1.5 M©, z = 0.00, = 
2.0 kms-S i^C/i^C = 10, e(C) = 8.24, e(N) = 8.35, e(0) = 
8.83. 



Fig. 1 reveals several discrepancies. As outlined in Sect. 
1, the origin of these discrepancies can be either a cali- 
bration problem or a problem related to the generation 
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of the synthetic spectrum. Precisely because this research 
involves both theoretical developments on the model spec- 
tra and calibration improvements of the spectral reduc- 
tion method, one has to be extremely careful not to con- 
fuse technical detector problems with astrophysical issues. 
Therefore, several points have to be taken into account 
when scrutinizing these discrepancies. 

First of all, one has to take into account that the cali- 
bration process is a delicate one. Due to the high flux level 
of the observations, dark-current subtraction does not in- 
fluence the final spectrum much, but memory effects affect 
the reliability of band-2 observations. Therefore, mainly 
the band-1 data were used for the determination of the 
stellar parameters. The final result is however very sensi- 
tive to the Relative Spectral Response Function (RSRF) 
used. The RSRF was intensively measured in the lab be- 
fore launch by use of a cryogenic blackbody source with 
temperatures between 30 and 300K. After the launch of 
ISO these intermediate RSRFs were corrected for broad- 
band residues detected by comparing reference spectral 
energy distributions (SEDs) — synthetic spectra (Kurucz 
1993; van der Bliek et al. 1998) and composites of M. Co- 
hen (Cohen et al. 1992, 1996a; Witteborn et al. 1999) — 
with SWS observations processed with these intermediate 
RSRFs using stellar candles of different IR spectral sig- 
nature. By using a cubic spline smoothing algorithm and 
sigmaclipping the risk was minimized of introducing false 
spectral features and one also avoids possible bias by spec- 
tral features of one type of source (Vandenbussche, 1999). 
The absolute flux level is determined by use of photomet- 
ric data and composite SEDs of various sources. The flux 
calibration accuracy is estimated to be 5 - 7% in band 1, 
7-12% in band 2a, 7-15% in band 2b and 11 - 25% 
in band 2c, the worst accuracy being located at the edges. 
The spectral cover from AO - M7 was highly necessary, 
since the different structures of the various selected stars 
make it possible to decide upon the quality of the used 
RSRFs. To test our findings, high-resolution observations 
of two stars were also studied in detail. For a discussion 
about the quality of the RSRFs, we refer to the forthcom- 
ing papers which will describe the other 15 stars in our 
sample. 

In order to elucidate problems with the theoretical at- 
mospheric structure, it is important to know the relative 
importance of the different molecules (see Fig. 4) and the 
influence of the stellar parameters on the total absorp- 
tion by these molecules (see Sect. 4). Fig. 4 shows clearly 
that especially the molecules CO, OH and SiO are help- 
ful in determining the atmospheric parameters. In spite 
of the moderate resolution of ISO-SWS, we will demon- 
strate that one can pin down the stellar parameters of 
cool giants very accurately from these data. This is due 
to the large wavelength range of ISO-SWS, where differ- 
ent important molecules determine the spectral signature. 
Since these different spectral features do each react in an- 
other way on a change of one of the several heterogeneous 



stellar parameters, it is possible to improve on the stellar 
parameters. Using a metallicity of z = —0.30 yields, e.g., 
a synthetic spectrum which has too low a flux density in 
the wavelength range from 2.6 - 2.8 /xm, with the line 
contrast I^/hv being too strong. Fig. 2 turned out to be 
extremely useful to ascertain the gravity: the inverse in- 
fluence of the gravity on the OH lines compared to the CO 
and SiO lines, together with the steep slope in band IB 
and band ID make it possible to pin down the logarithm 
of the gravity within 0.15 dex! Because the microturbu- 
lent velocity does not only act on the CO lines, but also 
on the strongest SiO and OH lines, one may disentangle a 
wrong microturbulence and a wrong carbon abundance. 

In order to achieve the highest possible agreement be- 
tween the data and the different synthetic spectra, a sta- 
tistical approach is proposed. A choice was made for the 
Kolmogorov-Smirnov test. This well-developed goodness- 
of-fit criterion is applicable for a broad range of compar- 
isons between two samples, where the kind of differences 
which occur between the samples can be very diverse. 
An elaborate discussion about the Kolmogorov-Smirnov 
statistics can be found in Pratt & Gibbons (1981) and 
Hajek (1969). In this goodness-of-fit test a parameter, Vfc, 
is defined in the present case as the ratio of the ISO-SWS 
flux to the synthetic flux at frequency point ■ These pa- 
rameters Vfe are assumed to be independent random vari- 
ables with E{Vk} — 1 and E{V^\ = u^, with a constant 
unspecified variance . The random variables 



n 

fe=i 



(i = 1, . . . ,n - 1) 



(3) 



(with n the total number of frequency points) are then 
asymptotically behaving like i/n and hence are asymptot- 
ically equidistributed on the interval [0, 1]. These variables 
Yi will be formally compared with the quantiles of the uni- 
form distribution F[x) = x, < a; < 1, using a modifica- 
tion of the Kolmogorov-Smirnov test. Taking the unspec- 
ified distribution function of the sample Yi, . . . , y„ as G, 
the null hypothesis is F = G. The Kolmogorov-Smirnov 
test is based on the maximum difference. Specifically, the 
test statistics are given by 



Am = max |G„(<) - -F„(t)|. 



(4) 



One rejects F = G if is 'too large' in comparison to 
an appropriate critical value. 

Practically, our modification involves the calculation 
of the value /3, defined as 



/3 = 



In sup 

l<i<n-l 



n 



(5) 



Firstly, we will compute the asymptotical distribution 
function of /3. By defining 

v.- 1 



X,, = 



(6) 
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one has that E{X,} = and E{Xf} = 1. So, 



P = v?i sup 

l<i<n-l 



n 



In sup 

l<i<n-l 



fc=i 



E 14 



(7) 



'n sup 

l<i<n-l 



(8) 



By defining k [nt] , where [nt] is the greatest integer less 
than or equal to nt, Eq. (8) becomes 



(5 = sup 

0<t<l 



fc=i 




fc=i 



77 ^ 



(9) 



Using the strong law of large numbers 
1 " 

lim — y Xk = 0, almost surely 



k=l 



and Donsker's theorem (Billingsley, 1968, p. 137) 

[nt] 



(10) 



(11) 



where W{t), < t < 1, is standard Brownian motion, i.e.. 
Wit) is a Gaussian process with W{Q) = 0, independent 
increments, and W{t) — W{s) has a normal distribution 
with mean zero and variance |i — s|, one sees that (3 con- 
verges in distribution to 



(3^ a sup \W{t)-tW{\)\. 

0<t<l 



(12) 



Here Wo{t) = W{t) - tW{l), < t < 1, is known as a 
Brownian Bridge. The definition of Wq implies that Wq is 
a Gaussian process with E{Wo{t)} = and 
i?{Wo(s)Wo(0} = ■s(l — t) if s < t. As a consequence 



r sup \W^{t)\, 

O' 0<t<l 



(13) 



which implies that the asymptotic distribution function of 
(5* is given by the Kolmogorov-Smirnov distribution 



P{(3* > A) = 2^(-l)*+icxp{-2i2A2} 



(14) 



(Billingsley, 1968, p. 105). For a test with significance level 
a one rejects the null hypothesis when (3* > \a, where 



2^(-l)'+iexp{-2*2Ai} = a. 



For a = 0.05 one has A^ = 1.36 and for a = 0.01 one 
has Aq, = 1.63. So, a good approximation to the level-a 
Kolmogorov-Smirnov test is to reject a fit if the empirical 
distribution function exits from the bounds 



y 



a;± A„(n- l)"^/^ < a; < 1 



(15) 



(see, e.g., Fig. 10.2 in BrockweU & Davis (1991)). 

The lower the /3-value, the better the accordance be- 
tween the observed data and synthetic spectrum. It is clear 
that (3 = a(3* gives us not only information about the cor- 
respondence between the observed and synthetic spectra, 
but also about a. Since T4 is the ratio of the ISO-SWS 
flux (Tohs) to the synthetic flux (Tsyn), it is obvious that 



and thus 



■p2 2 



E{J^ohs — -^^s; 
syn 



(16) 



(17) 



The unknown variance is thus an indication on the 
variance of the observed data. 

The Kolmogorov-Smirnov test globally checks the 
goodness of flt of the observed and synthetic spectra. An 
advantage of this test is that one very discrepant frequency 
point (e.g. due to a wrong oscillator strength) only mildly 
influences the final result. To avoid too high an influence 
of the less reliable calibration at the band edges, the pa- 
rameter f3 was computed for each sub-band separately. 
In each sub-band, several sub-intervals were considered 
to the main absorber in the relevant wavelength range, 
and the parameter /3 was computed for each of this sub- 
intervals. The resulting {Y — i^)-figures turned out to be 
very useful to reveal systematic problems (see, e.g.. Fig. 
5). 

Very low /?-values (< 0.09) were achieved, proving a 
good relation between observed and synthetic spectra. The 
objective /3-parameters, which are - in se - deviation esti- 
mating parameters, were then further used to improve the 
stellar parameters. A sensible rule of thumb would then be 
to select the synthetic spectrum with the lowest /3- values, 
among those that are consistent with the null hypothesis. 
In addition, it is a useful sensitivity analysis to compare 
various spectra with acceptable /3-values in key charac- 
teristics. The high flux accuracy enables us to distinguish 
between /3 = 0.15 and /3 = 0.03 (see Fig. 5), although both 
values are in the same region of accepting the null hypoth- 
esis. From the sensitivity analysis, it seemed to be useful to 
determine a maximum acceptable /3- value (i.e. a minimal 
deviation) as an objective criterion. A synthetic spectrum 
is accepted to represent the true stellar parameters of a 
Tau (and the other stars) when the /3-values are lower than 
the values given in column 2 of Table 3. These empirical 
/3max-values give, as already mentioned (Eq. (17)), an in- 
dication on the uncertainty on the observational spectrum 
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Fig. 5. Fig. (a-1) shows the comparison between band 1 of the ISO-SWS data of a Tau (black) and the synthetic spectrum 
(red) with stellar parameters Tetr = 3900 K, log g = 1.60, M = 1.5 Mq, z = 0.00, Ct = 2.0 kms-\ ^'^C/^^C = 10, e(C) = 8.24, 
e(N) = 8.35, e(0) = 8.83 and 6d = 20.67; in Fig. (a-2) the synthetic spectrum has as stellar parameters T^s = 3850 K, log g = 
1.50, M = 1.5 Mq, z = -0.15, = 2.0 kms^S ^^C/^^C = 10, e(C) = 8.24, e(N) = 8.35, e(0) = 8.83 and Od = 20.77 mas. In Fig. 
(b-1) and (b-2) the ratio of the ISO-SWS flux to the synthetic flux of respectively (a-1) and (a-2) is plotted. The difference 
between the partial sum Yi and the straight line i/n is plotted against the number i of the frequency point in Fig. (c-1) and 
Fig. (c-2). The parameter f3 of the different sub-bands is indicated in Fig. (b-1) and Fig. (b-2), whereas the parameter (3 of 
the entire band 1 is indicated in Fig. (c-1) and Fig. (c-2). 



determined by e.g. the pointing accuracy, the pointing jit- 
ter of ISO, uncertainties on the RSRFs, problematic dark- 
current subtraction, ... In Fig. 6 the standard deviation, 
obtained when rebinning the oversampled spectrum of a 
Tau to the expected resolution, of each bin is shown. Ob- 
viously, the ISO-SWS spectrum is less accurate in band 
lA than in the other bands. This is reflected in the higher 



maximum acceptable /3-value in Table 3. In column 3 of 
Table 3, the final /3-values for a Tau are given. 

Due to the smaller weight which are given automati- 
cally to small features, the traditional comparison between 
observed and synthetic spectra by eye-ball fitting is still 
necessary as a complement to this Kolmogorov-Smirnov 
method in order to reveal systematic errors in those fea- 
tures. The final error bars on the atmospheric parameters 
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Table 3. /3-value for the different sub-bands and molecules in 
that sub-band. In the second column the maximum acceptable 
value is given; in the third column the correspondent value for 
Q Tau. 



sub-band/ 


maximum acceptable 


/3- value 


molecule 


/3- value 


for a Tau 


la 


0.06 


0.040 


lb 


0.05 


0.034 


Id 


0.04 


0.036 


le 


0.04 


0.027 


CO (la) 


0.06 


0.040 


CO (lb) 


0.05 


0.039 


OH (lb) 


0.04 


0.024 


OH (Id) 


0.04 


0.036 


OH (le) 


0.03 


0.025 


SiO (le) 


0.03 


0.022 



STANDARD DEVIATION OF ALPHA TAU 




Fig. 6. Standard deviation of the different bins, obtained when 
rebinning the oversampled spectrum of a Tau to the expected 
resolution. 



tity Q depends on other independent parameters, e.g. a 
and b, the uncertainty s is estimated from the expressions 

or, equivalently 

The power of the Kohnogorov-Smirnov test has al- 
ready been proven in the field of astronomy, even multi- 
dimensional versions of the Kohirogorov-Smirnov test are 
developed for astronomical purposes (Fasano & France- 
schini 1987; Gosset 1987). Fasano & Franceschini (1987) 
demonstrate even that the Kolmogorov-Smirnov test is a 
much more goodness-of-fit test than the often used x^- 
tcst. It is, however, the first time that this test is used to 
evaluate synthetic spectra with respect to observational 
data. 

To summarize the strategy: 

1) from the large set of stars and two high-resolution ob- 
servations, it is possible to indicate calibration problems 
(sec forthcoming paper); 

2) in the comparison between observed and synthetic 
spectrum, some discrepancies were clearly visible. The 
knowledge on the relative importance of the different 
molecules (Fig. 4) and on the influence of the various 
stellar parameters (Fig. 2) makes it possible to elicit the 
origin of these discrepancies. This is due to the presence of 
different molecules in this large wavelength range, whose 
spectral shape and strength are a very good diagnostic 
tool to reveal the stellar parameters; 

3) once a high level of accuracy is achieved, the Kolmo- 
gorov-Smirnov test can be used for further refinement. 
From the {Y — F) plots, one can deduce systematic errors, 
yielding an indication on inaccurate stellar parameters. 
By using the results illustrated in Fig. 2, one can improve 
the stellar parameters. The stellar parameters, used to 
generate a synthetic spectrum, are taken to represent the 
'true' stellar parameters, once the /3- values are lower than 
a certain maximum value. 



are then estimated from 1. the intrinsic uncertainty on the 
synthetic spectrum (i.e. the possibihty to distinguish dif- 
ferent synthetic spectra at a specific resolution i.e. there 
should be a relevant difference in /3-values) which is thus 
dependent on both the resolving power of the observa- 
tion and the specific values of the fundamental parame- 
ters, 2. the uncertainty on the ISO-SWS spectrum which 
is directly related to the quality of the ISO-SWS observa- 
tion , 3. the value of the /3-parameters in the Kolmogorov- 
Smirnov test and 4. the still remaining discrepancies be- 
tween observed and synthetic spectrum. However, no ex- 
act formula can be given to compute the error bars due to 
the dependence of the several parameters. When a quan- 



6. Discussion 

Our best fit involves a synthetic spectrum with stellar pa- 
rameters Teff' = 3850 K, log g = 1.50, M = 2.3 Mq, z 
= -0.15, 6 = 1-7 kms-i, 12C/13C = 10, e(C) = 8.35, 
£(N) = 8.35, e(0) = 8.83 and 6*^ = 20.77 mas (see Fig. 
7). The uncertainties are estimated to be ATcfi = 70 
K, A log g= 0.15, A^t = 0.3 kms-i, Az ^ 0.20, Ae(C) 
= 0.20, A£(N) = 0.25, Ae(0) = 0.15, A^^C/i^C = 1 and 
A9d ~ 0.03. The quoted accuracy of 0.03 mas does not 
take the absolute photometric flux accuracy into account, 
which is 5 % in band 1, resulting in Ado < 0.83 mas. 
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Fig. 7. Comparison between band 1 and band 2 of the ISO-SWS data of a Tau (blaclc) and the synthetic spectrum (red) with 
stellar parameters Toff = 3850 K, log g = 1.50, M = 2.3 M©, z = -0.15, it = 1.7 l?ms-\ ^^C/"C = 10, e(C) = 8.35, e(N) = 
8.20, £(0) = 8.83 and Qd = 20.77 mas. The composite SED of Cohen et al. (1992) is plotted in blue. 



Due to the minor influence of the stellar mass, the un- 
certainty on this parameter, which follows from the con- 
struction of a synthetic spectrum, is quite large. But by 
using the parallax this uncertainty is constrained to AlVI 
= 0.8 Mq. The parallax measurement by Hipparcos for a 
Tau is 50.09±0.95 mas, which corresponds to a distance of 
19.96±0.39 pc. With an angular diameter 9d of 20.77±0.83 
mas, a radius R = 44.6 ±2.0 R0 is found, which combined 
with log g = 1.50±0.15, implies a mass of 2.3±0.8 Mq. A 
refinement in the mass determination was proposed by El 
Eid (1994), who noted a correlation between the ^^O/^^O 
ratio and stellar mass for evolved stars. Using this ratio in 
conjunction with evolutionary calculations. El Eid (1994) 
found a mass of 1.5 Mq for a Tau. A comparison of the 
deduced luminosity and effective temperature with evolu- 
tionary tracks (cfr. Mc William 1990) indicates a mass of 
2.0 Mq for a Tau. 

The accuracy of the method can be tested by making 
a grid around the resultant final stellar parameters. This 
is illustrated for a Tau in Table 4. The /3-values, corre- 
sponding to the final stellar paramters, for the different 
sub-bands are /3ia = 0.040, /3ib = 0.034, /3id = 0.036, 
Pie = 0.027 (see Table 3). In Table 4, the /3- values are 



given for a change of one of the main stellar parameters, 
being the temperature, the gravity, the metallicity, the 
carbon abundance, the nitrogen abundance and the oxy- 
gen abundance. From this table, one can deduce that there 
are four other acceptable sets of parameters, indicated by 
a grey box. The final stellar parameters — as mentioned 
above — were chosen to represent the 'true' atmospheric 
parameters of a Tau, due to the lower /3-value in band lA, 
where the important molecule CO is absorbing. 

According to Fig. 20 in Besscll et al. (1998), a Tau, 
with a (y — K) colour index of 3.67, has a bolometric cor- 
rection BCk of 2.60 ± 0.05. From the bolometric magni- 
tude TOboi = -0.21±0.05, the distance d = 19.964±0.379, 
the assumption that the interstellar extinction may be ne- 
glected, and by adapting the absolute bolometric magni- 
tude of the Sun Mboi,© = 4.74 (according to Bessell et al. 
1998), one deduces the luminosity of a Tau L = 382 ± 23 
L0. A radius R of 44.6 ± 2.4 Kq (calculated from the 
distance d and the angular diameter Od) yields an effec- 
tive temperature of 3828 ± 180 K. This is in very good 
agreement with the effective temperature deduced from 
the ISO-SWS spectrum, being Toff = 3850 ± 70 K. This 
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Table 4. /3-values from the Kolmogorov-Smirnov test for the 
evaluation of several synthetic spectra w.r.t. the ISO-SWS 
spectrum of a Tau. If not specified, the used stellar param- 
eters were Toft = 3850 K, log g = 1.50, = 1.7 kms"\ [Fe/H] 
= -0.15, e(C) = 8.35, e(N) = 8.35, e(0) = 8.83, 6^ = 20.77. 
The /3- values for these stellar parameters are given as /Sfinai- 









/3lA 


/3lB 


/3lD 


/3lE 






0.060 


0.050 


0.040 


0.040 




/3final 


0.040 


0.034 


0.036 


0.027 


Tcfi 




3710 


0.080 


0.050 


0.059 


0.027 






3780 


0.054 


0.028 


0.026 


0.017 






3920 


0.031 


0.043 


0.084 


0.028 






3990 


0.044 


0.046 


0.104 


0.041 


log g 






0.037 


0.093 


0.071 


0.015 






1.35 


0.034 


0.063 


0.064 


0.018 






1.65 


0.045 


0.028 


0.044 


0.027 






1.80 


0.051 


0.055 


0.032 


0.027 


z = 






0.037 


0.094 


0.075 


0.013 






0.30 


0.034 


0.062 


0.066 


0.015 






0.00 


0.045 


0.028 


0.043 


0.024 






0.15 


0.050 


0.054 


0.029 


0.030 


e{G) 




8.05 


0.047 


0.132 


0.040 


0.018 






8.20 


0.042 


0.065 


0.045 


0.021 






8.50 


0.036 


0.093 


0.063 


0.031 






8.65 


0.065 


0.174 


0.085 


0.044 


e(N) 




7.85 


0.043 


0.032 


0.051 


0.019 






8.10 


0.041 


0.031 


0.051 


0.021 






8.60 


0.034 


0.032 


0.055 


0.025 






8.85 


0.031 


0.034 


0.059 


0.028 


£(0) 




8.63 


0.057 


0.049 


0.096 


0.046 






8.73 


0.035 


0.040 


0.075 


0.034 






8.93 


0.052 


0.019 


0.026 


0.018 






8.03 


0.070 


0.043 


0.028 


0.019 



latter Tes, in conjunction with L = 382Lo, would give R 
= 44.05 ± 5.55 Rq and Qd = 20.53 ± 2.6 mas. 

It is clear that a gravity of log g = 1.50 ± 0.15 dex 
agrees better with a 'physical' gravity (log g=: 1.44 from M 
= 2 M0 and R = 44.6 R0) determined by the mass M, the 
radius R — or the effective temperature Teff and the lumi- 
nosity L — and/or evolutionary tracks than with a 'spec- 
troscopic' gravity of e.g. Lambert & Ries (1981), Kovacs 
(1983) and Luck & Challener (1995) (see Appendix A). 
It is well known that the accuracy of spectroscopic deter- 
minations of gravities from ionization equilibria or molec- 
ular equilibria for individual stars is not very good for 
red giants (cf.. e.g.. Trimble & Bell 1981, Brown ct al. 
1983, Smith & Lambert 1985). Smith & Lambert (1985) 



suggested that the low spectroscopic gravity log g = 0.8 
may be due to non-LTE effects on iron lines. After correc- 
tions for these effects. Smith & Lambert (1990) finally de- 
termined values for the various atmospheric parameters, 
with which our results agree very well. Our method for 
gravity determination from infrared spectra appears to be 
much less affected by NLTE effects, being based on a broad 
range of diagnostics over a large spectral interval. A great 
advantage of the used Kolmogorov-Smirnov test, is also 
that several deviation estimating parameters /3 are calcu- 
lated, which are not independent and are each of them 
related to the specific pattern of behaviour of absorption 
of the different molecules. E.g. the abundance of carbon is 
mainly determined from the band-lA and band-IB data, 
while a wrong gravity is traced from the /3iB,iD-values. 
Simulations with other (log g, £(C))-combinations never 
reached the same level of agreement as for the (log g = 
1.50, e(C) = 8.35)-combination. Using e.g. log g = 1.35 
requires a lower carbon abundance (than e(C) = 8.35) in 
order to fit the CO line strengths. The best fit then is ob- 
tained for £(C) = 8.25, with /3ia = 0.033 and /3ib = 0.024. 
The corresponding /3iD-value is then, however, 0.056 aris- 
ing from too steep a continuum. 

The mass M, luminosity L and radius R are the three 
fundamental parameters of a (spherical) star where two of 
these may be replaced by the surface gravity g = GM/R^ 
and the effective temperature (7Teff''=L/(47rR^). Both the 
mass and the luminosity are measurable physical quanti- 
ties, whereas the radius of a star is a fictitious quantity be- 
cause a star is a gaseous sphere and does not have a sharp 
edge. The relevant observable quantity is the center-to- 
limb variation (=clv) of intensity or limb darkening. The 
sun is the only star whose civ one can observe directly. Di- 
rect diameter determinations are based on a standard pro- 
cedure which evaluates interferometric, lunar occultation 
or binary eclipse data by determining a radius position on 
the base of a parameterized approximation (e.g. uniform 
disk) of a model-predicted limb-darkening curve (see e.g. 
Scholz 1998). The radius is then taken at /i = (/i cos 0, 
with 9 the angle between the radius vector and the line 
of sight). Published limb-darkening coefficients (for an 
overview, see Scholz (1998)) are then needed to convert 
to the limb-darkened corrected angular diameter. In Ta- 
ble 2, this limb-darkened angular diameter is listed for 
the direct angular diameter determinations, except when 
otherwise mentioned. The angular diameter value derived 
by us from the ISO-SWS data is however a 'theoretical' 
radius deduced from an indirect method. It corresponds 
to the definition crTcff^ = L/ (47rR(Ti.oss = 1)^) and log 
g = GM/R(rross = 1)^ (Plez et al. 1992). Baschek et al. 
(1991) have summarized various radius definitions. One of 
the most common theoretical definitions for the radius is 
the optical depth radius. Assuming that most of the fiux 
emerges from a layer with optical depth r w 1, is the base 
for defining a radius R^ = t^{tx ~ 1), which depends of 
course on the extinction coefficient at this wavelength. In 
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Fig. 2 of Scholz (1998) one can see clearly that there is no 
trivial correlation between the civ shape and the position 
of this ta = 1 radius. The angular diameter (and radius) 
obtained from the ISO-SWS data should therefore be only 
compared to other indirect spectrophotometric angular di- 
ameter values, resulting from e.g. comparing the observed 
flux at a distance d from the star with the surface flux, the 
surface brightness method (introduced by Barnes & Evans 
1976) or the infrared or Rayleigh- Jeans flux method (from 
Blackwell & Shallis 1977). The quadratic dependence of 
the observed flux on the angular diameter makes an accu- 
rate determination of 9d from the ISO-SWS data possible. 
With the exception of a few references given in the cata- 
logue of Fracassini et al. (1988), the angular diameter 9d 
= 20.77 ± 1.07 mas is in very good agreement with the 
other indirect determined angular diameter values. 

The carbon abundance is set by the strength of the 
CO absorption, while the strength of the OH lines yields 
the oxygen abundance. Due to the weakness of the NH 
and CN lines in a Tau (see Fig. 4), it was impossible to 
determine £(N) directly from the SWS spectrum. Since 
the chemical equilibrium detemines however the shape of 
the synthetic spectrum, a too inaccurate nitrogen abun- 
dance would result in high /3- values, which did not happen 
when using e(N) = 8.35 (see Table 4). The large depen- 
dence of the CO line- width on £(C) implies an uncertainty 
Ae(C) < 0.10 dex. This can also be seen in Table 4: chang- 
ing the carbon abundance by 0.15 dex results in e.g. too 
high a /3iB-value so that only for e(C) no other acceptable 
set of /3-values appears in Table 4. Due to, however, the 
larger standard deviation after rebinning in band lA and 
the beginning of band IB than in the other sub-bands, the 
total uncertainty on e(C) is taken to be 0.20 dex. From 
Table 4, one also can deduce that the error bars on e(N) 
and £(0), being 0.25 dex and 0.15 dex respectively, are 
realistic values. The carbon, nitrogen and oxygen abun- 
dance, together with the low isotopic ^^C/^^C ratio of a 
Tau are an argument for a Tau being on the first giant 
branch. 

With a luminosity of 382 Lg and Teft = 3850 K, a Tau 
lies slightly above the position where core helium-burning 
starts for stars with 1.0 Mq < M < 2.0 Mq and has an 
age of ~ 5 Gyr in the evolutionary diagram of Claret & 
Gimenez (1995). Consequently there are two possibilities 
for the evolutionary stage of a Tau: 

(i) a Tau is still on the first giant branch, or 

(ii) a Tau has ignited helium-burning and either had not 
enough time to move into the region of the clump stars or 
is already on the asymptotic giant branch. 

Following the discussion of Kovacs (1983), it seems im- 
possible to decide whether a Tau is in phase (i) or in 
phase (ii). So far, the low ^^C/^^C (= 10) value cannot be 
explained by standard evolutionary models for low-mass 
stars. Charbonnel (1994) introduced an extra-mixing pro- 
cess that occurs later on the giant branch and produces 
an additional decrease of the surface ^^C/^'^C and C/N 



values. This extra-mixing process should be only efficient 
when the hydrogen-burning shell has reached the discon- 
tinuity in molecular weight. Charbonnel (1994) held the 
process of Zahn (1992) — who proposed a consistent pic- 
ture of the interaction between meridional circulation and 
turbulence induced by rotation in stars — responsible for 
the required extra-mixing on the red giant branch of low 
mass stars. 

We want to stress that the current stellar parameters 
are obtained by using the 1-dimensional models as de- 
scribed in Sect. 4.1. Since there are indications of problems 
with the temperature distribution in the outermost layers 
of the theoretical models (see forthcoming paper), a fu- 
ture project is to extend the theoretical models to include 
the ability to conduct temperature perturbations in order 
to simulate convection, the presence of a chromosphere 
or a dust shell and a change in opacity. Using these im- 
proved models, it will be interesting to study the changes 
in the abundance (especially carbon, nitrogen and oxygen) 
needed to match the line strengths of CO (fundamental 
and first overtone)lines, SiO (fundamental and first over- 
tone) lines, H2O fundamental lines and OH fundamental 
lines in the ISO-SWS spectra and high-resolution Fourier 
Transform Spectrometer (FTS) spectra. In addition, by 
studying the effect of changing the microturbulcnce ^t, it 
might be possible to identify a typical region associated 
with a specific spectral type in the 3-dimensional (AT, 
^t, e)-plane, with e being the abundance of carbon, nitro- 
gen of oxygen. So far, it was still not possible to identify 
such a region from the many different values for the stel- 
lar parameters of red giants derived from different meth- 
ods and/or data (cf., e.g., Lambert & Ries (1981) versus 
McWilliam (1990)). Only a systematic study may reveal 
such a region. 

The composite SED of Cohen et al. (1992) is also plot- 
ted on Fig. 7 and hes ^ 2 % higher than the SWS data, 
which is well within the quoted flux calibration accuracy. 
Cohen et al. (1992) noted however a global accuracy of 
~ 1.5% for this a Tau SED and a local accuracy which 
is somewhat lower. The accuracy quoted by Cohen et al. 
(1992) thus appears slightly optimistic. 

The remaining discrepancies seen in band 2 are due 
to the memory effects which affect the reliability of ob- 
servations in this band. The obtained synthetic spectrum 
will be used as a test vehicle for the memory correction 
procedures which are currently under development. 

7. Conclusions 

The full scan AOTOl speed-4 ISO-SWS observation of 
a Tau has been used to determine the stellar parame- 
ters of this K5 giant. Due to the complementarity of the 
way parameters affect the spectrum in the wavelength 
range from 2.38 - 12 /im — like the presence of differ- 
ent molecules which are each dependent in another way 
on the several stellar parameters — it was possible to 
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pin down the atmospheric parameters with a high ac- 
curacy. This resulted in the following parameters Toff = 
3850 ±70 K, log g = 1.50 ±0.15 , M = 2.3 ±0.8 M©, z = 
-0.15±0.20, 6 = 1.7±0.3 kms-i, ^"^C/^^C = 10±1, £(C) 
= 8.35±0.20, £(N) = 8.35±0.25, e(0) = 8.83±0.15 and Od 
= 20. 77 ±0.83 mas. These parameters were compared with 
and discussed w.r.t. the results listed by other authors, one 
of the most striking result being that one should rely more 
upon a physically determined gravity than upon a spec- 
troscopic gravity. The good consistency between the stel- 
lar parameters deduced from the ISO-SWS data, colours, 
high-resolution (optical) spectra, IRFM, ... proves that it 
is possible to extract good quantative informations from 
these intermediate-resolution ISO-SWS spectra! 

Appendix A: Comments on the difTerent stellar 
parameters of a Tau. 

1. Mc William (1990) based his results on high-resolution 
spectroscopic observations with resolving power 40000. 
The effective temperature was determined by empirical 
and semi-empirical results found in the literature and from 
broad-band Johnson colours. The gravity was ascertained 
by using the well-known relation between g, Teff , the mass 
M and the luminosity L, where the mass was determined 
by locating the stars on theoretical evolutionary tracks. 
So, the computed gravity is fairly insensitive to errors in 
the adopted L. High-excitation iron lines were used for the 
metallicity [Fe/H] in order that the results are less spoiled 
by non-LTE effects. The author refrained from determin- 
ing the gravity in a spectroscopic way (i.e. by requiring 
that the abundance of neutral and ionized species yield the 
same abundance) because 'a gravity adopted by demand- 
ing that neutral and ionized lines give the same abundance, 
is known to yield temperatures which are ~ 200 K higher 
than found by other methods. This difference is thought to 
be due to non-LTE effects in Fe I lines. ' By requiring that 
the derived iron abundances, relative to the standard 72 
Cyg, were independent of equivalent width, the microtur- 
bulent velocity was found. 

2. Lambert & Ries (1981) have used high-resolution 
low-noise spectra. The parameters were ascertained by de- 
manding that the spectroscopic requirements (ionization 
balance, independence of the abundance of an ion versus 
the excitation potential and equivalent width) should be 
fulfilled. The effective temperature was found from the 
Fe I excitation temperature and the model atmosphere 
calibration of the excitation temperature as a function of 
Tefi. As quoted by Ries (1981), Harris et al. (1988) and 
Luck & Challener (1995) their Teft and log g are too high 
and should be lowered by 240 K and 0.40 dex, respec- 
tively. The isotopic ratio ^'^C/^^C was taken from Tomkin 
et al. (1975), while the luminosity was estimated from the 
K-line visual magnitude Mv{K) given by Wilson (1976) 
and the bolomctric correction BC by Gustafsson & Bell 
(1979). The abundances of carbon, nitrogen and oxygen 



were based on C2, [O I] and the red system CN lines re- 
spectively. Luck & Challener (1995) wondered whether the 
nitrogen abundance [N/Fe] = —0.20 quoted reflects a ty- 
pographic error and should rather be [N/Fe]=±0.20, re- 
sulting in £(N)=7.86, which is more in agreement with 
being a red giant branch star. 

3. The effective temperature Teff and the angular di- 
ameter 9d given by Blackwell et al. (1990) were determined 
by the infrared flux method (IRFM), a semi-empirical 
method which relies upon a theoretical calibration of in- 
frared bolometric corrections with effective temperature. 
One expects that the IRFM should yield results better 
than 1 % for the effective temperature and 2 - 3 % for 
the angular diameter. The final effective temperature is 
a weighted mean of T(J„), T(K„) and T(L„), with J„ at 
1.2467 nm, K„ at 2.2135 fim and L„ at 3.7825 fim. 

4. - 5. The stellar parameters quoted by Tsuji (1986, 
1991) were based on the results of Tsuji (1981), in which 
the temperature was determined by the IRFM method. 
A mass of 3 Mq was assumed to ascertain the gravity. 
Tsuji (1986) has used high-resolution FTS spectra of the 
CO (first overtone) lines to determine e(C) and (,t by as- 
suming that the abundance should be independent of the 
equivalent width of the lines. In Tsuji (1991) CO lines of 
the second overtone were used, where a standard analysis 
yielded results of e(C) = 8.39 and a linear analysis of weak 
lines led to £(C) = 8.31. 

6. High- resolution spectra of OH {Ai/ = 2), CO(Ai^ = 
3), CO{Av = 2) and CN{Av = 2 : A^U - A^S) were 
obtained by Smith & Lambert (1985). They have used 
{V — K) colours and the calibration provided by Ridgway 
et al. (1980) to determine Teff. Using the spectroscopic 
requirement that £(Fe I) = £(Fe II) yields log g = 0.8 dex, 
which is too low for a K5III giant. They suggested that 
the reason for this low value is the overionization of iron 
relative to the LTE situation. They then have computed 
the surface gravity by using a mass estimated from evo- 
lutionary tracks in the H-R diagram. The metallicity was 
taken from Kovacs (1983) and for the microturbulence 
they used Fe I, Ni I and Ti I lines, demanding that the 
abundances are independent of the equivalent width. Us- 
ing the molecular fines, they determined £(C), £(N), £(0) 
and 12C/13C. 

7. Harris & Lambert (1984) have taken Teff, log g and 
determined by Dominy, Hinkle & Lambert in 1984, a 

reference which we could not trace back. The isotopic ra- 
tio 12C/13C was adopted from Tomkin et al. (1975). The 
carbon abundance was found by fitting weak ^^C^^O lines 
at 1.6, 2.3 and 5 /.im. 

8. Kovacs (1983) obtained observations at the 1.52 m 
telescope of ESO at La Silla. By using different IR colour 
indices (y-R,V-I,V~J,V-K,V-L,R~MoT the 
broad-band photometry and 52 — 56, 52 — 63, 52 — 72, 52 — 
80, 52 — 86 for the narrow-band photometry) the effective 
temperature was found. The gravity and microturbulent 
velocity were determined in a spectroscopic way, where the 
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abundance derived from strong and medium lines should 
equal the abundance derived from weak lines for the right 
value of £,f The ^^C/^'^C ratio was taken from Tomkin et 
al. (1976). Using a parallax of 0."050 and 9d of 24 mas, 
resulted in a radius of 49 ± 4 Rq , which corresponds to 
a luminosity L of 463 ± 59 L0, while My = —0.7 and a 
bolometric correction taken from Johnson (1966) yields L 
= 413 Lq. 

9. Lambert et al. (1980) took model parameters from 
published papers which constrained the ^^C/^'^C ratio. 
The parameters for a Tau were based on Tomkin et al. 
(1976) and Lambert (1976). Tomkin et al. (1976) have 
used IR colours for the determination of Tcff and the mi- 
croturbulence was ascertained by fitting the theoretical 
curve of growth to the ^^CN curve of growth (^^CN(2-0) 
around 8000 A, weak ^'^CN{A-0) Hues around 6300A and 
weak i2CN(4-2) lines around 8430 A). Together with ^^CN 
lines around 8OOOA, the isotopic ratio ^^C/^^C computed. 
The gravities were estimated from the effective tempera- 
ture, the mass and the luminosity. 

10. By fitting the Engelke function (Engelke 1992) to 
the template of a Tau, Cohen et al. (1996b) could deter- 
mine Teff and 6d. A gravity of 2.0 was adopted, though in 
Cohen et al. (1992) a value of 1.5 — taken from Smith & 
Lambert (1985) — was used. 

11. By taking the mean value of the different temper- 
atures derived by calibrations with photometric indices 
{U - B, B V ~ R; V - I, V ~ J, V - K, V - L) 
Fernandez- Villacanas et al. (1990) have fixed the effective 
temperature. For the gravity, the DDO photometry in- 
dices C(45 - 48) and C(42 - 45) were used. The Fe I lines 
served for the determination of S,t ■ 

12. van Paradijs & Meurs (1974) have adopted the an- 
gular diameter from Currie et al. (1974). To determine the 
effective temperature they used several continuum data, 
the curve of growth with the line strengths of Fe I lines 
and the surface brightness. Requiring that the neutral and 
ionized lines of Fe, Cr, V, Ti and Sc gave the same abun- 
dance, yielded the gravity. Gravity, parallax and angular 
diameter yielded the mass, while the luminosity was de- 
termined from the effective temperature, the parallax and 
the angular diameter, van Paradijs & Meurs (1974) quoted 
that Wilson (1972) found [Fe/H] = -0.69. 

13. Tomkin & Lambert (1974) obtained photoelectric 
scans of the red CN lines with a resolution of 0.05 A for 
the red scans and 0.09 A for the infrared scans. The effec- 
tive temperature and gravity were adopted from Conti et 
al. (1967), who have ascertained the effective temperature 
from the absolute magnitude derived from the K emission- 
line width and the gravity from the effective temperature, 
the mass and the luminosity. The microturbulcnt velocity 
was determined using a curve of growth technique and the 
I2(^^i3(^ Y'dtio was obtained directly from the horizontal 
shifts between the curves of growth for ^^CN and ^'^CN 
lines. 



14. Also Luck & Challener (1995) have determined 
the effective temperature using photometric data (DDO 
C(42-45) and C(45-48), Geneva B2-V1, Johnson V-K, 
J — K, V — R and B — V). Two methods were used to fix 
the surface gravity. The first one determined the 'physi- 
cal gravity' by using the mass M and the radius R, with 
the mass M determined by Tos, L{Mv{K), EC) and the- 
oretical evolutionary tracks and the radius R by T^g and 
L{Mv{K), BC). The 'spectroscopic gravity' is based on 
the ionization balance of Fe I and Fe II lines. The dif- 
ference between these two gravities was very large, being 
1.0 dex! When the Fe II oscillator strengths were modi- 
fied to reflect a solar Fe abundance of 7.50 (instead of the 
used 7.67), then the spectroscopic gravity scale would rise 
by +0.25 dex (still far less than 1.0 dex) resulting in an 
increase of [Fe/H] of H-0.15 dex. They have quoted both 
pluses and minuses for the spectroscopic gravity and a 
single plus for the physical gravity. The microturbulence 
was ascertained by forcing no dependence of abundance 
(derived from individual Fe I lines) upon the equivalent 
width. In determining the carbon abundance using the C2 
Swan system lines and [C I] lines, the nitrogen abundance 
and isotopic ratio ^^C/^^C using CN lines and the oxygen 
abundance from [O I] lines, they always found a better 
trend with temperature when the spectroscopic gravity 
was used and also a better agreement between the two 
carbon indicators. When comparing their results with the 
ones of Lambert & Ries (1981), they found that a Tau was 
always the most discrepant star. Their very low spectro- 
scopic gravity, compared with the results of other authors 
using different methods and taking into account possible 
non-LTE effects, may be an indication that a value of log 
g 1.5 is in better agreement with the real gravity of a 
Tau. 

15. Aoki & Tsuji (1997) took the same stellar parame- 
ters as Tsuji (1986, 1991), but they now used CN-lines to 
determine e(N) and £,f 

16. Bonnell & Bell (1993) constructed a grid based on 
the effective temperature Tcs of Manduca et al. (1981). 
They used ground-based high-resolution FTS spectra of 
OH and [O I] lines. The requirement that the oxygen abun- 
dances determined from the [O I] and OH line widths agree 
amounts to finding the intersection of the loci of points de- 
fined by the measured widths in the ([0/H], log g) plane. 
For a Tau they found large discrepancies in the [O I] oxy- 
gen abundance, which leads to a spread of +0.3 dex in 
log g. The determination of was based on the OH-lines, 
but was hampered by a lack of weak lines from this radi- 
cal. Using a least square fit, they found = 1-05 kms^^ 
for the OH(Aj/ = 2) lines and = 1.62 kms'^ for the 
OH(Aj/ = 1) fines, refiecting that the OH(A^ = 2) se- 
quence lines are formed at greater average depth relative 
to the OH(Aj/ = 1) lines. For the [O I] lines a microtur- 
bulcnt velocity of 1.5 kms~^ or 2.0 kms~^ was used. 
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17. Ridgway et al. (1982) determined the limb-darke- 
ning-corrected angular diameter 0ld = Od using the lunar 
occultation technique. 

18a-b. Di Benedetto & Rabbia(1987) used Michelson 
interferometry by the two-telescope baseline located at 
CERGA. Combining this angular diameter with the bolo- 
metric flux Fboi (resulting from a directed integration us- 
ing the trapezoidal rule over the flux distribution curves, 
after taking interstellar absorption into account) they 
found an effective temperature of 3970 K, which is in good 
agreement with the results obtained from the lunar oc- 
cultation technique. Di Benedetto (1998) calibrated the 
surface brightness-colour correlation using a set of high- 
precision angular diameters measured by modern inter- 
ferometric techniques. The stellar sizes predicted by this 
correlation were then combined with bolometric flux mea- 
surements, in order to determine one-dimensional (T, V- 
K) temperature scales of dwarfs and giants. 

19. Mozurkewich et al. (1991) used the Marklll Optical 
Interferometer. The uniform-disk angular diameter ^ud 
had a residual of 1 % for the 800 nm observations and less 
than 3 % for the 450 nm observations. The limb-darkened 
diameter was then obtained by multiplying the uniform- 
disk angular diameter with a correction factor (using the 
quadratic limb-darkening coefficient from Manduca 1979). 

20. Quirrenbach et al. (1993) have determined the u- 
niform-disk angular diameter in the strong TiO band at 
712 nm and in a continuum band at 754 nm with the 
Marklll stellar interferometer on Mount Wilson. Because 
limb darkening is expected to be substantially larger in 
the visible than in the infrared, the measured uniform- 
disk diameters should be larger in the visible than in the 
infrared. This seems, however, not always to be the case. 
Using the same factor as Mozurkewich et al. (1991) we 
have converted their continuum uniform-disk value (19.80 
mas) into a limb-darkened angular diameter, yielding a 
value of 22.73 mas, with a systematic uncertainty in the 
limb-darkened angular diameter of the order of 1 % in ad- 
dition to the measurement uncertainty of the uniform-disk 
angular diameter (Davis 1997). 

21. Volk & Cohen (1989) mentioned a distance of 
20.0 ± 2.0 pc. The effective temperature was directly de- 
termined from the literature values of angular diameter 
measurements and total flux observations (also from lit- 
erature). The distance was taken from the Catalogue of 
Nearby Stars (Gliese, 1969) or from the Bright Star Cat- 
alogue (Hoffleit 1982). 

22. Blackwell et al. (1991) is a revision of Blackwell et 
al. (1990) where the H~ opacity has been improved. They 
investigated the effect of the improved H~ opacity on the 
IRFM temperature scale and derived angular diameters. 
Also here, the mean temperature is a weighted mean of the 
temperatures for J„, Kn and L„. Relative to Blackwell et 
al. (1990) there was a change of temperature up to 1.4% 
and an decrease by 3.5 % in Od- 



23. By using a TeS-{V — I) transformation of John- 
son (1966), Linsky & Ayres (1978) have determined the 
effective temperature. 

24. Bell & Gustafsson (1989) first determined the tem- 
perature from the Johnson K band at 2.2 /xm by use of 
the IRFM method. By comparison with temperatures de- 
duced from the colours Glass J— iJ , H — K, K — L and K] 
Cohen, Frogcl and Pcrsson J — H,H — K,K — L and K; 
Johnson V-J, V-K, V-L and A'; Cousins V~R, R-I; 
Johnson and Mitchell 13-colour and Wing's near-infrared 
eight-colour photometry, they found that Teff(IRFM) was 
~ 80 K too high, by which they corrected the tempera- 
ture. The gravity and the metallicity were adopted from 
Kovacs (1983). 

25. Taylor (1999) prepared a catalogue of tempera- 
tures and [Fe/H] averages for evolved G and K giants. 
This catalogue is available at CDS via anonymous ftp to 
cdsarc.u-strasbg.fr 

26. Burnashev (1983) has determined Tgff, log g and 
[Fe/H] from narrow-band photometric colours in the visi- 
ble part of the electromagnetic spectrum. 

27. Fracassini et al. (1988) have made a catalogue of 
stellar apparent diameters and/or absolute radii, listing 
12055 diameters for 7255 stars. Only the most extreme 
values are listed. References and remarks to the different 
values of the angular diameter and radius may be found in 
this catalogue. Also here these angular diameter values are 
given in italic mode when determined from direct meth- 
ods and in normal mode for indirect (spectrophotometric) 
determinations. 

28. Perrin et al. (1998) have derived the effective tem- 
peratures for nine giant stars from diameter determina- 
tions at 2.2 /xm with the FLUOR beam combiner on the 
IOTA interferometer. This yielded the uniform-disk an- 
gular diameter of a Boo and a Tau of our sample. The 
averaging effect of a uniform model leads to an underesti- 
mation of the diameter of the star. Therefore, they have 
fitted their data with limb-darkened disk models published 
in the literature. The average result is a ratio between the 
uniform and the limb-darkened disk diameters of 1.035 
with a dispersion of 0.01. This ratio could then also be 
used for the uniform-disk angular diameters of 7 Dra and 
(3 Peg, fisted by Di Benedetto & Rabbia (1987), and a Get, 
which was based on a photometric estimate. Several pho- 
tometric sources were used to determine the bolometric 
flux, which then, in conjunction with the limb-darkened 
diameter, yielded the effective temperature. 

29. Engelke (1992) has derived a two-parameter ana- 
lytical expression approximating the long-wavelength (2 
~ 60 /xm) infrared continuum of stellar calibration stan- 
dards. This generalized result is written in the form of a 
Planck function with a brightness temperature that is a 
function of both observing wavelength and effective tem- 
perature. This function is then fitted to the best empirical 
flux data available, providing thus the effective tempera- 
ture and the angular diameter. 
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30. Blackwell & Shallis (1977) have described the In- 
frared Flux Method (IRFM) to determine the stellar an- 
gular diameters and effective temperatures from absolute 
infrared photometry. For 28 stars (including a Car, a Boo, 
a Cma, a Lyr, /? Peg, a Cen A, a Tau and 7 Dra) the an- 
gular diameters are deduced. Only for the first four stars 
the corresponding effective temperatures are computed. 

31. Scargle & Strecker (1979) have compared the ob- 
served infrared flux curves of cool stars with theoreti- 
cal predictions in order to assess the model atmospheres 
and to derive useful stellar parameters. This comparison 
yielded the effective temperature (determined from flux 
curve shape alone) and the angular diameter (determined 
from the magnitudes of the fluxes). The overall uncer- 
tainty in Tcft is probably about 150 K, which translates 
into about a 9 % error in the angular diameter. 

32. Manduca et al. (1981) have compared absolute flux 
measurements in the 2.5 - 5.5 yum region with fluxes com- 
puted for model stellar atmospheres. The stellar angular 
diameters obtained from fitting the fluxes at 3.5 /zm are 
in good agreement with observational values and with an- 
gular diameters deduced from the relation between visual 
surface brightness and {V — R) colour. The temperatures 
obtained from the shape of the flux curves are in satisfac- 
tory agreement with other temperature estimates. Since 
the average error is expected to be well within 10%, the 
error for the angular diameter is estimated to be in the 
order of 5 %. 

33. Smith & Lambert (1990) have determined the 
chemical composition of a sample of M, MS and S giants. 
The use of a slightly different set of lines for the molecular 
vibra- rotational lines, along with improved ^/-values for 
CN and NH, lead to a small difference in the carbon, nitro- 
gen and oxygen abundance compared to Smith & Lambert 
(1985). 

Appendix B: The approximate absorption 
coefficient TZ 

Kjaergaard et al. (1982) examined two cases for which they 
derived equations for the approximate absorption coeffi- 
cient TZ = l,y/Ku. In case (a) most of the carbon and nitro- 
gen are still considered to be free atoms, assumed for Tcff 
> 5000 K, while in case (b) almost all the C and N are in 
the form of CO and N2, assumed for Tcff < 4500 K. Look- 
ing at the relatively high partial pressure of CN, it seems 
that the assumption of complete association of N into N2 
is not valid for our giant models. Attributing all depletion 
to the formation of N2 seems to be more valid in a dwarf 
model (Bell & Tripicco 1991). The formula can however 
still give a qualitative view of the influence of the differ- 
ent parameters. For case (b) Kjasrgaard et al. deduced an 
approach for CN: 



7^(CN) 



Using the approach of Kjaergaard et al., similar equations 
for CO, SiO, OH, H2O and NH are derived. 
For CO, we have 



7^(C0) 



^ k(CO) 
" «(Hi) 



n(CO) 

_ n(C l)n(0 I)K(H-) 
~ K(CO)n(H I)n(c) ' 

In the assumption that all carbon is locked into CO, 
7?.(C0) may also be written as 



(B.l) 
(B.2) 
(B.3) 



7^(C0) cx 



n(C) 



n(H-) 
_ n(C)K(H-) 
n(H I)n(e) 
« (C/H)(Fc/H 



-0.4„-0.5 



(B.4) 

(B.5) 
(B.6) 



with the approximation of n(c) and n(H)/n(e) in the line- 
forming region being (Kjasrgaard et al., 1982) 



n(e) Z (Fe/H)°■^g°■^ and 
n(H)/n(e) °S {Fe/R)-^. 



(B.7) 
(B.8) 



The approximate absorption coefficient of SiO is de- 
fined as 



7^(SiO) 



^ K(SiO) 
" «(Hi) 



cx 



n(SiO) 
n(H-) 

n(Si l)n(0 I)K(H-) 



(B.9) 
(B.IO) 
(B.ll) 



K(SiO)n(H I)n(c) 

For oxygen-rich giants, n(0 I) may be approximated by 

n(0 I) w n(0) - n(CO) w n(0) - n(C). 

The equilibrium of silicon is dominated by SiO and Si I 
(Si II is present in the deeper layers as well). For example, 
in the model with Teff= 4050 K, log g = 1.00, z = 0.00, 
M = 1.5 Mq, = 2.0 kms"\ the proportion of the 
number densities are given in Table B.l. Since the ratio 
n(SiO)/n(Si I) changes rapidly in the outer layers of the 
photosphere, a rough approximation to n(Si I), valid in 
the region where —2 < log(ro) < —1, is then 

n(Si I) « (Si/H)n(H I), 
so that 



n{SiO) cx (Si/H)(0/H - C/H 



(N/H)"-"(C/H)(0/H - C/H)-^(Fe/H)-"-^g 



-0.1„.-0.75 



n(H I) K(H-) 
n(e) K(SiO) 
cx (Si/H)(0/H - C/H)(Fe/H)-i. 



(B.12) 
(B.13) 
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Table B.l. Proportion of number densities for silicon for the 
model with parameters Teff — 4050 K and log g = 1.00. 



log(ro) 


SiO 


Si I 


Si II 


-5 


83 


17 





-4 


54 


46 





-3 


26 


74 





-2 


9 


90 


1 


-1 


3 


93 


4 








35 


65 



7^(0H) 



Using the same notations, TZ{OH) is written as 
k(OH) 



n(OH) 
oc — -, 

n(H-) 

_ n(0 I)n(H I)K(H-) 
~ K(OH)n(H I)n(e) 

mm I)K(H-) 

cx (0/H-C/H)(Fc/H)-i. 
Also for H2O we obtain 

-<«^«>^^ 

n(H20) 
°^ n(H-) 

_ n(H2)n(0 I)K(H-) 
~ K(H20)n(H I)n(e) 
_ n(H I)n(H l)n(0 I)K(H-) 
" K(H2)K(H20)n(H I)n(e) 

„,0/H-C/H,^m 

= (0/H - C/H)(Fe/H)-i-6g° ^ 

The approximate absorption coefficient of NH is 
ten as 



7^(NH) 



k(NH) 



oc 



n(NH) 

n(N I)n(H I)K(H-) 
K(NH)n(H I)n(e) 



B.14) 

B.15) 

B.16) 

B.17) 
B.18) 

B.19) 

B.20) 

B.21) 

B.22) 

B.23) 
B.24) 
writ- 

B.25) 
B.26) 
B.27) 



The equilibrium of the nitrogen species is dominated by 
the N2 formation. So, 

n(N 1)2 = n(N2)K(N2), 



with 

n(N2) « 0.5(N/H)n(H I), 
resulting in 

n(N I) « [(N/H)n(H I)K(N2)]^/^ 
Consequently, 

7^(NH) cx (N/H)°-5n(H I)"-^n(e)-^ 
0.5 n(H 1)0-5 



cx (N/H) 



-n(e) 



-0.5 



n(e)0-5 

oc (N/H)°-5(Fe/H)-°-5(Fe/H)-''-2g 
= (N/H)0-5(Fe/H)-°-7g-o-25. 



-0.25 



(B.28) 

(B.29) 

(B.30) 
(B.31) 
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